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Abstract 


Pharmacophore-based virtual screening is an effective, inexpensive and fast approach to discovering useful starting points for drug dis- 
covery. In this study, we developed a pharmacophore model for the main proteinase of severe acute respiratory syndrome coronavirus (SARS- 
CoV). Then we used this pharmacophore model to search NCI 3D database including 250, 251 compounds and identified 30 existing drugs 
containing the pharmacophore query. Among them are six compounds that already exhibited anti-SARS-CoV activity experimentally. This 
means that our pharmacophore model can lead to the discovery of potent anti-SARS-CoV inhibitors or promising lead compounds for further 


SARS-CoV main proteinase inhibitor development. 
© 2004 Elsevier SAS. All rights reserved. 
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1. Introduction 


The infection of the newly emerged severe acute respira- 
tory syndrome coronavirus (SARS-CoV) is characterized by 
acute flu-like symptoms that progress to acute lung injury or 
acute respiratory distress syndrome with over 10% of mor- 
tality [1]. To date there are no universally recommended 
therapy for the disease. Many scientists are now making 
efforts to develop effective drugs against SARS. The combi- 
nation therapy of corticosteroid with lopinavir, ribavirin and 
ritonavir can improve clinical response and reduce mortality 
rates apparently [2,3]. Cinatl et al. [4] found that ribavirin, 
azauridine, pyrazofurin and glycyrrhizin are active against 
SARS-CoV. Barnard et al. [5] reported that calpain inhibitors 
and B-D-N*-hydroxycytidine exhibit inhibitory effects on 
SARS-CoV. 

Structure-based drug design focuses on two important 
approaches: one is receptor-based docking technique, another 
is pharmacophore-based virtual screening technique. A phar- 
macophore is the 3D arrangement of atoms or functional 
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groups essential for the compound to bind to a specific recep- 
tor [6]. The power of a pharmacophore model is to discover 
new leads by using 3D database pharmacophore searching 
and guide chemists to synthesize new compounds [7]. Sucha 
pharmacophore-based method has been successfully applied 
to many drug development programs [8—18]. The main pro- 
teinase of SARS-CoV plays an important role in virus repli- 
cation and is the primary target for drugs. The aim of this 
study is to develop 3D pharmacophore models for SARS- 
CoV main proteinase and expect to provide useful knowl- 
edge for anti-SARS drug design. 


2. Material and methods 


There are two methods to derive a reasonable pharmacoph- 
ore model. One is from the crystal structure of protein— 
ligand complex, another is based on molecular modeling of 
enzyme with its potential inhibitors. Here we used the experi- 
mental structure of SARS-CoV main proteinase complexed 
with its peptide inhibitor CMK (PDB ID 1UK4) [19] and the 
predicted structures of SARS-CoV main proteinase with six 
drugs/compounds [20] for establishing pharmacophore mod- 
els. The structures of CMK peptide and six compounds are 
shown in Fig. 1. 
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Fig. 1. The peptide and compounds used for pharmacophore generation. 


The POCKET module in LigBuilder program [21] was 
employed to obtain the pharmacophore models of SARS- 
CoV main proteinase. This approach was successfully ap- 
plied to the identification of novel inhibitors for alanine 
racemase [22]. The proposed pharmacophore model is a 
binding-site-derived pharmacophore model, which includes 
the following pharmacophore features of ligands binding to 
the enzyme’s active site: a positively charged nitrogen atom 
(ammonium cation) to represent a hydrogen bond donor 
(HBD), a negatively charged oxygen atom (as in a carboxyl 
group) to represent a hydrogen bond acceptor (HBA), and a 
carbon atom (methane) to represent a hydrophobic center 
(HPC). A pharmacophore model is generated for each pro- 
tein—ligand complex. 


3. Results and discussion 


Using CMK peptide and six compounds in Fig. 1, we 
generated a set of seven eight-point pharmacophore models 
for SARS-CoV main proteinase, which is listed in Table 1. 
These hypotheses exhibit different features due to the diver- 
sities of the compounds involved. Based on these models, we 
extracted a common four-point pharmacophore distance pat- 
tern shown in Fig. 2, Pl is HBA, HBD and HPC, P2 is HBA 


and HPC, P3 is HBA and HBD, P41s HBA and HBD. Sucha 
pharmacophore distance pattern was subsequently used for 
the 3D database search. 

The pharmacophore searching in 3D database was con- 
ducted in the 3D NCI database, which has 250, 251 open 
structures ready for searching. Table 2 summarizes the re- 
sults for similarity search of above-mentioned four-point 
pharmacophore model with constraints: (1) the distance 
ranges for P1P2, P1P3, P2P3, PI1P4, P2P4 and P3P4 (Fig. 2) 
are 2—3, 2—3, 1.5—2.5, 5—6, 3-4 and 3-4 A, respectively; (2) 
the compound is drug; (3) the antiviral probability is over 
70%; (4) the compounds labeled as “No Name” is not in- 


Table 1 
The eight-point pharmacophore models obtained by LigBuilder from seven 
peptide and compounds 


Pharma- Peptide/ Features 


cophore compounds 


CMK peptide HBD HBD HPC HBD HPC HBD HPC HBA 
Lopinavir HBA HBD HBA HBD HBA HBA HBA HBA 
Ritonavir HPC HPC HBD HBD HBD HBA HBD HBA 


Niclosamide HBA HBD HBA HBA HBA HBD HPC HBA 


YANN BP WN 


Promazine HBA HBD HBA HPC HBA HBA HBA HBA 
PNU HBA HPC HBA HBD HBD HBA HBA HBA 
UC2 HPC HBD HBD HBA HBD HPC HPC HBA 


HBA = hydrogen bond acceptor, HBD = hydrogen bond donor, HPC = hy- 
drophobic center. 


X.W. Zhang et al. / European Journal of Medicinal Chemistry 40 (2005) 57-62 59 


5.695 


3.534 


P4 


Fig. 2. Four-point pharmacophore distance pattern for SARS-CoV main 
proteinase. Here, Pl is HBA, HBD and HPC, P2 is HBA and HPC, P3 is 
HBA and HBD, P4 is HBA and HBD. 


Table 2 
Summary of NCI database search by four-point pharmacophores 
Pharmacophore features in Fig. 2 Hits 

Pi P2 P3 P4 
HBD HBA HBD HBD 987 
HBA HBA HBA HBA 794 
HPC HBA HBA HBD 286 
HPC HPC HBA HBD 305 
HPC HBA HBD HBD 298 
HBA HBA HBD HBD 881 
HBA HBA HBA HBD 895 


HBA = hydrogen bond acceptor, HBD = hydrogen bond donor, HPC = hy- 
drophobic center. 


cluded. After review of these hitlists, 30 drugs were selected 
for further analysis, their chemical structures and bioactivi- 
ties (documented in anti-HIV and anti-opportunistic infec- 
tion chemical compound database, which contains approxi- 
mately 100,000 compounds, http://www.apps1.niaid.nih. 
gov/struct_search/an/an_search.htm) are shown in Fig. 3 and 
Table 3, respectively. It is noted that almost all drugs exhib- 
ited anti-HIV activity, and some of them have activities 
against Mycobacterium tuberculosis (puromycin), influenza 


Bioactivity documented in HIV/OI therapeutics database 
(http://www.apps 1 .niaid.nih.gov/struct_search/an/an_search.htm) 

HIV 

HIV 

HIV, HSV, human cytomegalovirus, vaccinia virus, Cowpox virus 

HIV, HSV, human cytomegalovirus, varicella-zoster virus, vaccinia virus, 
COWpox virus 


HIV, HSV, vaccinia virus 

HIV 

HIV, human cytomegalovirus, HSV 
HIV 

HIV, M. tuberculosis 

HIV 

Human cytomegalovirus 

HIV, HSV, varicella-zoster virus 


HIV 
HIV, vaccinia virus, West Nile virus 
HIV, varicella-zoster virus 


HIV, HSV, human cytomegalovirus, varicella-zoster virus, measles virus 
HIV, HSV, vaccinia virus 

HIV, cowpox virus, vaccinia virus 

HIV, HBV 

HIV, HBV 

HSV 

HIV, HSV, HCV, influenza virus, dengue virus, measles virus, respiratory 
syncytial virus, rhinovirus, polio virus, vaccinia virus, Cowpox virus 
West Nile virus, cowpox virus, vaccinia virus, dengue virus, Japanese 
encephalitis virus, Yellow fever virus 

Cowpox virus, vaccinia virus, HSV, varicella-zoster virus 

HIV, HSV, cowpox virus, vaccinia virus, varicella-zoster virus 

HIV, vaccinia virus, human cytomegalovirus 


Table 3 

Thirty drugs obtained by four-point pharmacophore search in NCI 3D database 
Name NSC number Formula 
PD-ADI 218321 C4 HieNsO, 
Coformycin 277817 C,,H1.N.0; 
Zidovudine 602670 CigH 3N.0, 
Vira-A 404241 Ci 9H,3N50, 
Angustmycin C 53104 C,,H,;N;0; 
ARA-AMP 259272 C,9H,4N;0,7P 
Cordycepin 63984 CioH3N503 
Triciribine 154020 Ci3Hi6N.Ou 
5-AZCR 102816 C,H,.N,0; 
Puromycin 3055 C,~H39N70; 
CHETOMIN 28949] C,,H39N,0,S.4 
Vengicide 99843 C,.H,,N;0, 
Spongothymidin 68929 C1 9H14N,0, 
Arauridine 68928 C,H,,N,0, 
P-Ara-C 135962 C,5H,3N3;0, 
Pyrazofurin 143095 CyH,3N30, 
Thymidin 21548 C,9H,4N,05 
Radibud 38297 CoH, ,BrN,O; 
Alexan 63878 CoH, 3N30; 
Floxuridin 27640 C,H, ,FN,O. 
Gemcitabine 613327 C,H,,CIF,N;0, 
Dideoxy guanosine 619072 C,9H13,N;03 
2’,3’-Dideoxycytidine 606170 CaH,.N,0; 
5-Bromo-2’-deoxycytidine 61765 CoH,,BrN.0, 
Ribavirin 163039 C,H,.N,05 
Azauridine 32074 C,H, N30, 
Fialuridine 678514 C,H, FIN,O; 
Emanil 39661 CoH, ,IN,O; 
Tubercidin 56408 C,,H,4N,0, 
Viroptic 75520 Ci0H11F,N.05 


HSV, cowpox virus, vaccinia virus 
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Fig. 3. Chemical structures of 30 drugs obtained by four-point pharmaco- 
phore search in NCI 3D database. 


virus (5-bromo-2’-deoxycytidine), hepatitis virus (dideox- 
yguanosine, 2’,3’-dideoxycytidine and ribavirin), dengue vi- 
rus (azauridine and ribavirin), rhinovirus and poliovirus (rib- 
avirin). This should be very meaningful in consideration of 


Fig. 4. The mappings of six compounds that experimentally exhibited anti- 
SARS-CoV activity into the four-point pharmacophore model: (A) azauri- 
dine, (B) 5-bromo-2’-deoxycytidine, (C) dideoxycytidine, (D) dideoxygua- 
nosine, (E) pyrazofurin, (F) ribavirin. 


the following facts: (1) there are some links between SARS- 
CoV and HIV and HBV [23-25]; (2) similar structure pat- 
terns exist in SARS-CoV main proteinase with rhinovirus 3c 
protease, poliovirus 3c proteinase, HAV 3c protease, HCV 
Ns3 protease and dengue virus Ns3 protease [26]; (3) SARS- 
CoV has clinically similar symptoms with influenza virus/M. 
tuberculosis, such as fever, cough, pains, pneumonia and 
death [27]. 

Indeed, among the 30 drugs are six compounds that al- 
ready exhibited anti-SARS-CoV activity experimentally: 
azauridine, pyrazofurin, ribavirin, 2’,3’-dideoxycytidine, 
dideoxyguanosine, and 5-bromo-2’-deoxycytidine [4,5]. 
This shows that our pharmacophore model can lead to the 
discovery of potent anti-SARS-CoV inhibitors or at least 
provide some useful clues. Fig. 4 shows the mappings of the 
six compounds into the four-point pharmacophore model, 
which are mapped to 1-3 HBA, 1—2 HBD and 0-1 HPC. In 
addition, most of the remaining compounds have remarkable 
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Fig. 5. Superpositions of anti-SARS-CoV compounds with other com- 
pounds. (A) azauridine with 5-AZCR, alexan, arauridine and spongothymi- 
din; (B) 5-bromo-2’-deoxycytidine with email, fialuridine, floxuridin, radi- 
bud, thymidin and viroptic; (C) dideoxyguanosine with angustmycin C, 
cordycepin, tubercidin, vengicide and vira-A. 


similarities with one of the above six compounds, for ex- 
ample, azauridine with 5-AZCR, alexan, arauridine and 
spongothymidin; 5-bromo-2’-deoxycytidine with emanil, fi- 
aluridine, floxuridin, radibud, thymidin and viroptic; dideox- 
yguanosine with angustmycin C, cordycepin, tubercidin, 
vengicide and vira-A. The superpositions for these com- 
pounds are shown in Fig. 5. In summary, our results indicate 
that the existing 30 drugs identified by our pharmacophore 
model could be potential inhibitors against SARS-CoV, or at 
least good lead compounds for anti-SARS-CoV drug design. 
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